A planar extrapolation of the correlation problem that permits pairing 
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It was observed previously that the SU{N) extension of the Hubbard model is dominated by 
planar diagrams at large A'', but the possibility of superconducting pairing got lost in this extrap- 
- - - olation. To allow for pairing, we replace SU{N) by U{N,q), the unitary group in a vector space 

, of quaternions. At the level of the free energy, the difference between the SU{N) and U{N,q) 

extrapolations appears only to first nonleading order in A''. 
PACS Numbers; 71.27.+a 71.10.-w 71.10.Fd 
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I. INTRODUCTION AND MOTIVATION 



, The Hubbard model y belongs to a class of models describing strongly correlated electrons in systems such as, for 
example, high Tc ceramics, transition metal oxides and organic conductors However, in spite of much progress 
made over the last ~ 40 years, it is still difficult to determine the low temperature properties of such model systems 
1^ ' as a function of their chemical input parameters. 

Recently it was noticed that an SU{N) deformation of such models is dominated, for ^ 1, by planar diagrams 
in the sense of 't Hooft [| and the sum of the leading planar diagrams turned out to be very similar to the " Fluctuation 
Exchange Approximation (FLEX)" pioneered by Bickers and Scalapino However, unlike the FLEX approximation, 
, the new method permits a systematic construction of the generating functional in powers of 

' To exploit this new perspective on the correlation problem, we should study the evolution of the properties of a given 
I model as one moves from iV = oo to iV = 2. But precisely such a study is impossible in the SU{N) extrapolation of 
' the Hubbard model since for TV > 2 it contains no analog of the superconducting order parameter < ^a{x)sa0ipi3{y) >■ 
^ I The very feature that constituted one of the motivations of studying Hubbard like models in the first place got lost 
in this extrapolation. 

Some time ago it was recognized that U {N, q), the unitary group acting in a vector space of quaternions, provides 
an interesting way of implementing the large N limit because it generalizes the antisymmetric tensor Sa[3 in a natural 
way. Because of U{1, q) = SU{2), an extrapolation via U{N, q) with N — 1,2,.. is just as good an extrapolation as 
SU{N) with = 2, 3.. Results of what probably amounted to an U{N, q) extrapolation of the tJ model were given 
recently 0], but without any calculational details. 
[ We will show below that there is a unique U{N,q) extrapolation of the Hubbard model that permits pairing and, 
T— I ■ simultaneously, a topological expansion in which the planar diagrams dominate for A^ ^ 1. 

o ■ 
o : 

O . II. U{N, Q) INVARIANT INTERACTIONS BETWEEN BOSONS AND FERMIONS 

» ' 

a • „ 

' It is reasonable to begin our discussion by recalling the definition of the group U{N, q) Q. It acts in a vector space 
^ , of quaternions: 

; 

C ! i'a= i'a,0 + i ^ 'lpa,r'y'^ and t/i* = -0a,o - i X! ^a.rtx'' (1) 

r=1..3 r=1..3 

Tpa^i real for a = 1..N, fi = 0..3, 7?" = Pauli matrices 



The ipa^i are real because we wish to think of the matrices icr'', r — 1..3 as square roots of —1. The transformation 
matrix g is also quaternion valued, g = g'^ + with g^, fi = 0..3 real in order for the image to be a quaternion. 

As a further constraint, it must leave a quaternion valued length invariant: 

E E {r9+)a{9^P)a-^g+g = l (2) 

a=l..N a=l..N 

The properties of the Pauli matrices if imply g'^ = (^g'^ — ilf'd'^ and g^ = g^ — ilt~^ from which we deduce 
g'^ = a'^g^a^. In conjunction with g^g = 1 this latter relation implies g^cr'^g — ■ The relations 

g^ g = 1 and g'^cr^g = cP' (3) 

can be restated by saying that U{N , q) has two invariant tensors, SabO'aP ^^'^ ^ab^ap- 
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The quaternions ipa {^'a)^^ of eq(|]) are not all independent because they can be parametrized by only 4 real 
quantities for each value of a. To get rid of the redundant quantities, while keeping the attractive transformation 
properties of the quaternions, we define "conventional" complex electron operators 



^aa= J2 i^a).rir ^ ^ = [I ] (4) 

T = l,2 



that have only one spinor index. 
A nonvanishing average such as 



< iJaa{x)al^My) >+ (5) 



breaks charge but conserves U {N, q) and pairing is therefore possible in an U{N, q) extension of the Hubbard model. 

Following the strategy described in we will introduce auxiliary boson fields via a Hubbard Stratonovich trans- 
formation. For this we need (i) an invariant quadratic form in the auxiliary bosons and (ii) an invariant interaction 
between the auxiliary bosons and some fermion bilinears. Both requirements are met if the auxiliary fields are infinites- 
imal generators of C/(iV, q), so we consider them now. Writing finite transformations as g = expiz , z = zzq + cr* 
with complex, for the moment, and using g+g ~ 1 and a^g = we find 

z+ = z and z^ a'^ +a'^z^O (6) 

The first of these equations implies {zq — —zq, — ^} while the second means |z(f = — -zq, = Both taken 

together imply that z^ is real so that iz is a quaternion valued matrix (and therefore g = expiz also a quaternion). 
The preceding argument shows that U {N, q) depends on a total of ^N{N — 1) -f ^N{N + 1) = 2N^ + N independent 
parameters. In particular, for A'^ = 1, there are 3 generators, the antisymmetric matrix zq is absent and from the 
explicit paramctrization g = exp i z' a' we conclude that U{1, q) ~ SU{2). This conclusion is also very natural in view 
of the definition of the group according to cq(|^). The quadratic form in the Lie algebra 

I |z| |2 = ^Tr{zz+) = ^Tr (*zo + z^Tf ) {~iz+ + ^+lf) (7) 



rr(zoz+ + z^ z*+) = Y ""ab^ 



is positive and invariant and an invariant coupling between fermions and bosons is given by (iz*^ + '^'a') "0. Because 
of its symmetry properties, z^ couples only to roughly half the bilinears in ■(/'"''and ip in this expression: 

V'+ {iZo + z'ct') Ip^Y ^°b^ [i^ta^ba 6)] + ^ ab^ [V'aa"^«/3V'fc/3 + {a ^ b)] = ^ Z^j, J^^ (8) 

with Jab = \ \^toS^afi^b!i + h.c. and = {i, ct"} 

Indeed, symmetrization with respect to Hermitian conjugation amounts to symmetrization with respect to the indices 
a, h. 

III. HAMILTONIAN AND PARTITION FUNCTION 

The infinitesimal generators z can be converted into dynamical bosons z(r) with the help of a Hubbard Stratonovich 
transformation: 

const = J Dz exp -jyrJ2\«b + Jab)^ (9) 
const ■ exp dr-^ JabJab = _/ ^ 51 2 + ^""abJab) 



It remains to find out what energy of interaction is represented by the expression — i ^ Jab^ab- This is done in eq(|2 
of the appendix and leads to 
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X! -^abJab + 7ab7afc = ^ X! {^l^a) {^b^b) " {i't ^^1) ii^bSipb) + quadratic terms (10) 

3n|n| + quadratic terms 

where ipbS'ipb means ^pba^iTTipbT and e — ia^. As extrapolation of the Hamiltonian we choose, therefore 



6N 



(11) 



This Hamiltonian collapses to the conventional expression for = 1 and, as we shall see below, it scales correctly in 
A'' and has a simple topological expansion. 

By appealing to the Trotter technique, the partition fmiction of the Hubbard model may then be represented in 
path integral form: 



Z = Trexp-l3{Ho + H,nt) 

= const J D^Dz exp - J dr ^P*{dt + ho)^ + ^ ^ab^Lb + ""^ab^Lb 



(12) 



Here /lo is the hopping matrix of the model and we have used the fact noted above that z^^^ couples only to the 
symmetrized part of the current J^^jf,- A more compact and esthetically pleasing form of the action is: 



^ = |^^V*(at + /io)V^ + ^5]rr(zz+) + y^^^+zV'j withz = ^a'^;: 



(13) 



IV. FEYNMAN DIAGRAMS AND 'T HOOFT'S TOPOLOGICAL EXPANSION 



From eq(|l2|) and using either textbook knowledge Q or auxiliary currents, we find the following bare propagators 
of the non interacting theory: 

< r{V'ao(l)V'?/3(2)} >0= 5c.p5ab {^Q^^ (1.2) ^"'^ P-'^^S^^^*™ 5^p5ab9{h2) (14) 

<T{zZbS^)zZbS''))>^=5^u,.5{l,2)Pi:^,^^^^^^ 

Kb^a.b. = \ [5a,.aJb,M + (" 1)'^'° Jfc,,,, ] , p^^ . p^^ = p^^ 

The projector P'^ reflects the difference in symmetry properties between zq (antisymmetric) and Zi (symmetric) for 
i = 1,2,3. Above we have also indicated that the tensor structure of the full propagators is the same as that of the 
bare propagators, because of symmetry considerations pO| ]. There is also a vertex 

^aa4a6<^a/3V'h/3 (15) 

where a particle hole excitation splits into its constituents. Propagators and vertex are graphically represented in 
figure 1 
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FIG. 1. Single line electrons, double line bosons and the vertex in the U{N,q) model 

New, relative to the corresponding SU (N) diagrams, is the presence of two terms in the boson propagator. Repre- 
senting fermion propagators by oriented hnes, boson double lines inherit an orientation from the electron lines that 
flow through a vertex. This inherited orientation, however, is not respected by the first term of the propagator in 
eq(Oh and a consequence of this is illustrated in figure 2. 




FIG. 2. The boson propagator spoils the loop orientation 

As in SU{N) Yang Mills theory |] or in the SU{N) Hubbard model y a diagram can be viewed as a collection 
of index loops that are identified across double lines representing bosonic propagators. Each index loop spans a two 
dimensional area and in this way a topological surface gets associated with each diagram. By counting the number 
of vertices and index summations from each loop, we find that each diagram carries a weight (in powers of N) of 

weight = N'oops-^ ^ ^s.-^ s^^J^^ 
X — S2 — Si + Sq 

where S2 is the number of loops or area like pieces, 5*1 the number of propagators and Sq the number of vertices and 
where Si = |S'o relates the number of propagators and vertices. In the present U{N, q) extrapolation, non orientable 
surfaces contribute and a more general expression for the Euler characteristic must be used: 

X = 2-2h-c-o 

Here h stands for handles, c for crosscaps and o for openings (see e.g. for a detailed discussion of the classification 
of surfaces). A "crosscap" is the presence, in a loop, of a pair of identified segments, that are traversed with the 
same orientation. For large N, diagrams of the disk topology h — = c, o=l dominate and these can be easily 
summed (the sums of the nonleading diagrams for x = and x = — 1 are more complicated and will not be given 
here). 

V. COUPLED INTEGRAL EQUATIONS FOR THE LEADING GRAPHS 

Let us compute the electronic self energy to second order in perturbation theory: 
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with Zaa\bi3 = •^^ab'^'Q/j- From eq(^6|) of the appendix we find 

SaiQi,b2/32 (Ij 2) = (5ai62<5ai/32^(l; 2) 

-d, 2) = ^^^^^(1, 2)5(1, 2) f 2),(1, 2) 

Strictly speaking, this is just lowest order perturbation theory. But to leading order in TV we get all the rainbow 
diagrams and these are summed by rendering the outermost rainbow explicit and we do indeed obtain the last equation 
with the full propagators (to leading order in TV), but without vertex corrections (that are nonleading in iV). We 
now compute the self energy of the bosons to leading order (our notation was indicated in eq(|l^): 

n:^j::fc,(l,2) = ^ < {^l.K[,,i^b,0.) (1) • (2) > (16) 

To lowest order this gives, upon using the invariant tensor properties of the fermion propagator 

n!:,Vla2\2(l'2) = 7r(l,2)(T'^^'^-)^^,^„^,^ (17) 

( ^^'"''^)aifcia2fc2 — ^bia2'5fc2ai'^MlM2(~l) 
III 

^1,2) = -— G(1,2)G(2,1) 

It is shown in eq ( |30| ) of the appendix that the conventional Dyson equation for the propagator leads to a corresponding 
scalar Dyson equation for d: 



The coupled integral equations 



d = ——^ and g = (18) 

, , (5$ , , , 2 (5$ 

^(1'2) = -TT^TTT ^(1,2) = 



5.g(2,l) ' ' ' NSd{2,l) 

<i> = -|5(l,2)d(l, 2)5(2,1) 

must be solved self consistently. Using the differential form of Dyson's equations 



<7{l,2)^~^-^^Tr{logg-g/go) (19) 
7^(1, 2) = -^^^Tr (log d - d/do) 



and eq(|l|) we find that the functional 

N 

n = '^> + Tr (log g - 5/50) - (log d - d/do) 

is stationary with respect to variations in g and d. To understand the factor ^ j we note that the theory contains 2 A'^ 
complex fermions and ~ 2N^ real bosons or ^ times as many more complex bosons than fermions. 

This leading generating functional is of the same form as a corresponding leading functional of the SU (N) theory 
that was written down previously which is consistent with a general idea that different groups have the same 
leading asymptotics, once their dimensions are taken large. Differences appear, however, at the first nonleading order, 
because of the nonorientable part of the boson propagator. 
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VI. CONCLUSIONS 



The main results of this paper are the foUowing: 

• we constructed an J7(-/V, q) extrapolation of the correlation problem that is dominated by planar diagrams at 
large N and which permits pairing - 

• the surfaces associated with the diagrams of this U{N, q) theory now include non orientable ones - 

• at the leading planar level, the integral equations and generating functional are of the same form as in the 
SU{N) case. 

It is very encouraging that the planar diagram approach to the correlation problem permits pairing, once an 
appropriate extrapolation is chosen. This is absolutely essential for its future use in investigating low temperature 
superconducting order in models of doped ceramics, organic conductors and other systems of interest. By using 
quaternion like objects, one has kept one of the crucial features of the electrons, their spinor character. 

Regarding magnetic ordering, it is unclear at present whether antiferromagnetism at half filling will persist also for 
iV > 1 in this extrapolation. It is conceivable that the magnetization La^^)^, with La a basis in the space of 
generators of U{N,q), arranges itself to be antiparallel (in the sense of the metric gAB — TtLaL^) at neighboring 
points a; on a bipartite lattice, in a similar fashion as discussed in 

It remains to find out, whether the small parameter 1/N introduced here is "small enough" to confer predictive 
power to our approach. So far the only favorable indication we have is the overall agreement of the leading diagrams 
with those of FLEX plus the many interesting results obtained in this approximation . If the present approach is 
successful, it will sharpen the FLEX method and a number of tough problems in the domaine of strongly correlated 
electrons will become transparent. 
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VII. APPENDIX 



A. Calculation of Jl^^J^^^ 



We use the definition J^^ = i i^ta'^ap^bp + h.c. with cr^ = {z, ct*} and calculate in two steps via J^bJab = 

Jab^ab + 'J'ab^ab- Fust wc determine JabJab ■ 

E -^abJab = - i E [^ar^"- " ^ ^)] t^-^^-^ " ^ &)] (20) 

Second we determine ^ ~y ab~J ab ■ 

E '^ab^ab = E ap^bp [^t^'^^SlpbS + [a ^ b)] (21) 

= E ^^ia^bp [V'iyV'M + (a ^ b)] {2SasSpj - 5cp^^&) 

= E - (V^a ^fc) V-fc) - {^Ua) {^t-^b) - \ {^t^b) (V'a+^h) - \ {i'ti'b) i^P^^a) 



The symbol = means modulo quadratic terms that correspond to an extra chemical potential. Combining eqs( poy2l| ) 
we obtain 
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E -^abJab + ^ab^ab = ^ " i^ai'a) " 2 {^t^b) (V'.+ V'b) (22) 

Because there are only two spinor components we have {ipa'^b) (V'^V's) = ~ \ {'i't^'^t) ('0f)£V'&)- Using this on eq([2^) 
we find the required relation 

E ''abJlb + '^ab'Jab = E (^a eV-^) (V'ftEV'b) " (V'^V'q) {-^l^t '^b) + quadratic terms (23) 



B. Calculation of < T { [^(1)^(1)]^^^^ [^+(2)z(2)] J > 



It is convenient to organize the algebra as follows 

< t{[z(1)VX1)],^„^ [^+(2)z(2)]^^^J >=< T{ (1) (2)} ><T{V'(l)6,ftV'+(2)a2a2} > (24) 

< T{z(l)z(2)} > is linearly related to < T {z^^^^^ (l)z^^^(j^ (2)} > for which we have an expression in eqjl^) of the 
body of the paper: 

<T{z,,„,„,^,(l)z,,„,|,,^,(2)} >=Y.<P.<P.<T{<,bMXbM > (25) 

= d(l,2)(5ai,a2'5bi,62 • 2^ai/3l^a2/32 + <^i'^^'^)^o.^.b2K■a2 2 '^aift^a2/32 

with a^' = {i, o^}. We use the relations jcr^i^iCr^^^^ = Sa^a^ep^p^ and ^o-^j/JiCT^^/j, (-l)''^'" = '^ai/32'5/3ia2 in eqdH) 
to obtain: 

< 7'{^aiQi|6i/3i(l)^a2Q2|b2/32(2)} >—d{l,2) ((Sqi ,62 "^aift ' (^bi ,03 (S/3ia2 " (^ai ,02 £qiQ2 ' '^bl :f)2 ^ftft ) 

where we recognize the basic invariant tensors of the group U{N,q). Invoking the tensor character of the fermion 
propagator in eq(p^) , the expression in eq(p^ reduces to 

<T{za,c.^\b^f3^il)Za,c.,\b,f3A^)} >< T {^{1)^,0,^+ i2)a,a,} > (26) 

= rf(l,2)g(l,2) {2N + l)5a,M5c,p, 
which is the desired expression. 

C. Scalar form of the bosonic Dyson equation 

It remains to see how this self energy modifies the scalar coefficient -0(1, 2). We must evaluate the sum 

D = Do + DoWDo + ... = E (n^o)" = nn = ^-i^ tt (27) 

Df) is proportional to a projector P onto matrices of definite symmetry type, its inverse D^^ does not exist and, at 
first glance, Dyson's equation appears ill defined. In any case we need a Dyson type equation for the scalar quantity 
d that encapsulates the information of D. The quantities entering the series in eq(p7|) are 



i?(l,2) = 5^,,^,d(l,2)F^ "projection" (28) 

nMi,A'2 = 7r(l,2)(5^,^3(-l)'^^'°T^ "transposition" 



with 



P'' =\ + (-l)*"'"'5ai,62'^fcl 

= (-I)*^''4ia2'^b2ai "transposition- 
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Our definitions were chosen so that the matrices are contracted with each other in the same order as one writes them. 
The relations 

come as no surprise and they mean that the operator n^i-^^ g^^i^g g^g multiplication by the scalar 7r(l,2). We can 
now sum up the series eq(p7[) 

D = Do + DoIlDo + ... = Pdo V (Trdof = Pdo-^—r = —f— (30) 
where {irdoY^ meant in the sense of convolutions and the last equation is the desired scalar Dyson equation. 
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